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Abstract 



Wave energy harvesting could be a substantial renewable energy source without impact 
on the global climate and ecology, yet practical attempts have struggled with problems of 
wear and catastrophic failure. An innovative technology for ocean wave energy harvesting 
was recently proposed, based on the use of soft capacitors. This study presents a realis- 
tic theoretical and numerical model for the quantitative characterization of this harvesting 
method. Parameter regions with optimal behavior are found, and novel material descriptors 
are determined which simplify analysis dramatically. The characteristics of currently avail- 
able material are evaluated, and found to merit a very conservative estimate of 10 years for 
raw material cost recovery. 
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INTRODUCTION 



The problem of adequately supplying the world with clean, renewable energy is among 
the most urgent today. It is crucial to evaluate alternatives to conventional techniques. One 
possibility is energy harvesting from ocean waves, which has been proposed as a means of 
offsetting a large portion of the world's electrical energy demands |1| . However, the practical 
implementation of wave energy harvesting has met with obstacles, and the development 
of new methods is necessary [2j. Oceanic waves have large amplitude fluctuations that 
cause devices to fail due to excessive wear or during storms. A strategy to overcome these 
catastrophic events could be to base the harvesting mechanisms on soft materials. 

Soft, stretchable rubber capacitors are possible candidates for energy harvesting |3]- 
6j, that have already been tested in a realistic ocean setting [8]. They were originally 
introduced as actuators [9HT2], capable of high actuation strains of more than 100% and 
stresses of more than 1 MPa. With a soft capacitor, mechanical energy can be used to 
pump charges from a low electrical potential U to a higher, such that the electrical energy 
difference can be harvested [3j. This is made possible by the large changes of capacitance 
under mechanical deformation. Although the method is simple and proven [3HH], it is still 
not clear to which extent the approach is practically useful, which is the concern of this 
paper. Of the many electro-active polymers, it appears that soft capacitors could have the 
highest energy densities |13| . 

For the purpose of a broad and realistic investigation, a minimal yet realistic model is 
proposed that takes into account the mechanical and dielectric properties of the soft capacitor 
material, including losses and limiting criteria. The model also includes cyclic mechanical 
driving, and an electrical control mechanism consisting of a switchable electrical circuit. The 
quality of the energy harvesting is characterized by efficiency and gain measures, which were 
evaluated for simulations on a very large set of varied system and material parameters. 

The so far only commercially available soft capacitor material, DEAP (Danfoss PolyPower 
A/S), is used here as benchmark [HHT6]. DEAP consists of a sheet of silicone elastomer with 
smart compliant metal electrodes and is a realistic candidate with adequate properties for 
energy harvesting |17| . The high voltage required for operation of soft capacitor actuators 
has led to research efforts aiming at lowering the voltages by modifying their mechanical or 
dielectric properties. Preferably, they should have low mechanical stiffness and high dielectric 
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constant, e r [T81 - I22) . In general, the dielectric constant can be adjusted from 2 to more than 
1000, however, at higher values such strategies usually cause excessive conductivity and 
premature electrical breakdown. 



THE MODEL 

The model describes the dynamics of the deformations of and the voltage across the soft 
capacitor. It includes mechanical and dielectric properties. Losses appear mainly electrically 
in the electrodes and in the charging and discharging circuits. The external force is sinu- 
soidal and linearly coupled to the length of the soft capacitor; we regard these as minimum 
requirements for the description of coupling to near-coast oceanic surface waves. Technically, 
elaborate wave-capacitor coupling mechanisms are possible, yet here the focus is placed on 
the general principle. Mechanical conversion efficiency and electrical energy gain factors are 
calculated from time-integrated losses, mechanical input and electrical output. 

The deformation of the polymer film is described in terms of the deformation ratios 
Xi = Li/ (with i = 1, 2, 3 for the ), where Li and L\ are the final and the initial 

dimensions, respectively. The electrical field is applied in z-direction and the mechanical 
driving force fx acts in x-direction. For elastomers, which are amorphous and non-crystalline, 
wolume conservation can be assumed, A1A2A3 = 1. Also, the pure shear condition is assumed, 
A2 = constant. 

The time-dependent mechanical response is chosen as a Kelvin-like fluid [23j, 




which balances the internal material stress Tu, the hydrostatic pressure p, the external 
forces fi, and the viscous damping characterized by 7. Eliminating the pressure from both 
equations, utilizing the volume constraint (A3 = A^f 1 ), and splitting the material stress into 
a mechanical and an electrical part yields 




Tu=p + - 7Aj 



(1) 
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For o"Mech; the Yeoh model [25] is chosen, with parameters Ciq, C20, C30: 
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FIG. 1: The electrical circuits. Left: charging process, right: discharging process. 

where A = \\ + A^ 2 — 2. The Young's modulus of the Yeoh material under pure shear 
conditions is Y = 8C\o. 

The electrical part of the internal stress 0"Eiec> also known as Maxwell stress is OEiec — 
—E r EoE 2 , where E = U/L 3 is the electrical field between the electrodes [U |2H 26J. It 
describes the stress due to the voltage difference U between the capacitor electrodes. The 
permittivity of an elastomer does not change at the relatively low levels of electrical field 
encountered here, and it also does not change when it is deformed, as has been repeatedly 
verified experimentally by researchers [24, 27J. 

The dynamics of the voltage between the electrodes can be derived from Q = CU, where 
Q is its charge and C the capacitance. Building the time derivative and rearranging this 
equation yields 



—U + -^ 1 -2U— 
Eqe t CAi \\ 

The first term corresponds to leakage current due to the finite polymer conductivity o"dc> the 
second accounts for the electrical control circuit, and the third for mechanical deformation. 
C = ErEoLiLt}/ L$f is the capacitance of the unperturbed capacitor, and C = C'\\ that of 
the deformed. 

Eqs. ([2]) and ^ define a system of two coupled ordinary differential equations (ODEs) 
for the stretch A and the voltage U. The only unknowns are the external driving force fx 
and the electrical current from the charging and discharging circuit. The force is chosen by 
a simple periodic driving /1 = /o + /a sin(o;t + if), where the actual geometry of the pick-up 
mechanism can help in adjusting bias force /o and amplitude f\. The angular frequency 
u = 27r/tc is determined by the period of the oscillating force, tc- 

A very simple charge and discharge control mechanism is chosen here. It is described by 
three points in time. The voltage is ramped on the soft capacitor at t = 0, stopped at t 1( 
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while discharging starts at t?,. The discharging process ends at tc- A phase difference, ip, 
adjusts the onset of the control relative to the driving force. 

The capacitor is charged during (i , ti) in an electrical circuit shown in Fig. [T] Kirchhoff 's 
laws state that Ui n + U + Uce + Ub. Iu = and that Ic = const. U\ n is the driving voltage of the 
external voltage supply, Uqe, and Ur Iu are the voltage drops over the capacitor electrodes 
and the charging resistor, with resistance Rqe and Ri n . This yields the current 

-RcE + -Rin 

which follows the voltage ramp Uj n (t) = ^-Ui njiaax , where Ui a>max is the maximal charging 
voltage. 

In the interval (£1,^2) the circuit is detached. Mechanical deformation and variation in 
voltage takes place due to the external variation in force. The response is influenced by 
viscoelastic damping, while current leakage leads to partial charge loss. 

During (t2,tc) the discharging electrical circuit is attached, cf. Fig. [Tj Now, no external 
voltage is applied and the discharging resistor Ro ut is used: 

Ic = - E " • (6) 

-KCE + -KOut 

After discharging the cycle starts over. 

In the ODEs Q and Q four relaxation processes with different time scales are given: 
one mechanical, tm, two electrical, tq and td, and rpc describing the loss of charge through 
the polymer material. In addition, there is the period of the driving, tc- The mechanical 
relaxation of the polymer to an equilibrium state is described by r^ 1 = dF(Ai)/dAi|A*- 
F(X) is given by ^ via Ai = F(Ai). For the Yeoh model this time scale is approximately 
tm = 0.55 s. The electrical relaxation time scales during charging and discharging are 
Tc = (Rce + RinjC and td = (Rce + Rout)C- They are typically in the order of 10 -5 s. The 
largest time scale, which we term the Maxwell time, describes the loss of charge through the 
polymer material rpc = ^j-. 

Material failure sets limits which are monitored during simulations and which have been 
taken into account: i) Ai < 5 avoids rupture, ii) T n > ensures a taut material, iii) 
E = < 20 ^ avoids intrinsic electrical breakdown, and iv) ^§fp-Af < H(\i) avoids the 
the electromechanical instability, analogous to pull-in failure [28, 29J. 

The mechanical work is generally defined as W = L[ J f(Xi,t)Xidt, where f(X,t) is an 
arbitrary force acting on the film. Using this equation and writing all terms in Eq. ^ as 
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TABLE I: Default parameters for energy harvesting simulations. 



Material and capacitor parameters 

Li, L' 2 , L' 3 0.1m, 0.2m, 5- 10~ 5 m 

Cio, C 20 , C* 30 139840 Pa, -6570 Pa, 1057 Pa 

e r 3.1 

io- 13 4 

r C e iooo n 

7 10 6 s/Pa 
Driving parameters 

A = 1N, /p = 2N, t c = 10s 

Circuit parameters 

^In,Max 651V 

Rin, Rom W 6 Q, 3.83 • 10 7 ft 

*i, *2, *c* = 27r/a; 
Limiting parameters 

•^Max 5 

E Max 20 Vm- 1 



forces allows the calculation of the work, for example WMech = L[ f Ai<j Mcch \idt. Energy 
is conserved over one cycle, hence WVisco + ^Driving = W^Mech + ^Maxwell- The electrical 
work is similarly W = fUIdt; again, several contributions should be evaluated. Energy 
conservation leads to Wcapacitor + Wu,in + W^R.in + W^ce = for the charging, and to 
^Capacitor + W^R.Out + W^r.ce = for the discharging. It is emphasized that the definition of 
the voltage allows for negative or positive values for the work done. 

The quality of the harvesting process is evaluated by the harvesting efficiency rj and the 
gain G. The efficiency rj is the ratio of electrical output energy to the total mechanical and 
electrical input energies (note that W^in.Eicc < by convention): 

^Driving — H^U.In 
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FIG. 2: An example of a cycle with the default values listed in Table [IJ The phase of the driving 
force is ip = 1.75929 and the switching times are t\ = 1.8 s and t2 = 4.32 s. 

The gain G is the net relative electrical energy gained compared to electrical energy invested: 

W R n ut 1 



G 



U.In 



(8) 



G is positive (negative) if electrical energy is gained (lost). Losses occur due viscoelasticity, 
leakage currents and finite resistance of capacitor electrodes and charging circuits. Both 
measures are essential to feasibility evaluation of energy conversion in practice. 

Eqs. @ and Q can only be solved numerically. The ratio of smallest to largest time 
scales is on the order of 10 6 , making the system stiff and solvable by methods such as the 
Radau, Rosenbrock, or Runge-Kutta solvers [30| I3T1J . Here, the Runge-Kutta 4th order 
method is used, with an integration time step smaller than the smallest time scale of the 
system. 



SIMULATION, OPTIMIZATION AND RESULTS 

Now follows a presentation of results obtained for the variation of the parameters of the 
harvesting cycle, the electrical circuit, and the materials. Some such parameters could be 
varied easily in an operating wave power plant, while others would be fixed by technical and 
design choices. 

The voltage-stretch curve for a simulation using the parameters listed in Table [I] are 
shown in Fig. [2] This particular example corresponds to fully optimized parameters for 
the DEAP material, which were found through procedures described in the sections below. 
The mechanical and dielectric parameters are taken from measurements performed on the 
DEAP soft capacitor material, and are consistent with the corresponding data sheets [15J. 
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FIG. 3: The harvesting efficiency rj and the gain G in dependence of the switching times t\ and 
t%. The surfaces are triangular because t\ < t%. The phase is fixed (p = 1.75929. 

The phase of the driving is ip = 1.75929, the switching times are t\ = 1.8 s and t 2 = 4.32 s. 
The efficiency in this particular example is rj = 0.42 while the gain is G = 0.47. The total 
harvested energy is W^out = 5.5 mJ. Compared to the mass of the material (the volume is 
L^LgLg = 10~ 6 m 3 , while the density is 1100 kg m" 3 ), the specific harvested energy for one 
cycle becomes 5Jkg _1 . Of course, this is small compared to energy density of 1500 Jkg -1 
(from which applies to the special case of prestretched polyacrylate glue. However, the 
result obtained here is more realistic; it is found by considering a full cycle, maximizing r] 
and with a limitation of the electrical field to a reasonable level of 20 Vm" 1 . 

Cycle parameter optimization 

The harvesting cycle allows a variation of the switching times £ 1; ti and the phase (p. 
Their impact on the optimal efficiency 77 and gain G is shown in Fig. [3] for the DEAP 
material (see Tab. [I]); the phase is fixed in this plot to ip = 1.76. As expected, efficiency and 
gain depend strongly on the cycle parameters. Note also that regions of maximum gain and 
maximum efficiency do not overlap perfectly. Hence, an optimization strategy for efficiency 
could result in zero energy gain and vice versa. 

In the following, both efficiency and gain are optimized in dependence of the circuit 
and material parameters. For each parameter set, ti, t 2 and (p are varied to maximize the 
efficiency and the gain, either by a Monte Carlo method using random sampling or the 
simplex method |31| . 
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Electrical circuit optimization 

In a power plant setting, the parameters L^Max and -Rout are adjustable. So is Ri n , 
though this would typically increase loss, and thus is not studied further here. Fig. [4] 
illustrates the result of this optimization. It was obtained by first choosing random values for 
E' = t/jr^Max/Lg and -Rout in the ranges shown, then optimizing the efficiency by adjusting 
£1, t2, via the simplex method. If the breakdown criteria are violated, the particular 
data point is plotted in red. The material properties are again taken from Tab. [Tj Clearly, 
the efficiency increases if the driving voltage increases. It is bounded by the breakdown 
criteria, and it is relatively unaffected by the output resistance in the investigated range. 
The maximum observed harvesting efficiency was about 0.42. 



Dielectric material parameter optimization 

Fig. [5] shows the numerically optimized efficiency for varying polymer conductivity o~dc 
and permittivity e Te \ for two different values of the stiffness of the polymer. The stiffness 
is varied by scaling of C±o, C20, and C30. For each individual parameter set (e r , croc) 1000 
randomized sets of values of Uj n , Rout, t\, ti and ip were generated; other parameter are 
held fixed, cf. Table |lj Out of these data sets, the one with the highest value of 77 (or G) not 
violating the limit criteria was determined and plotted (therefore, the observed 77 and G do 
not correspond to identical Ui n , i?out, ti, t<i and ip values). 

The results in Fig. [5] confirm intuition: First, r) generally decreases with o"dc and increases 

0.75 1 
0.50 - 
0.25 - 

io 1 " Ls ' s ^5 

10 s 

FIG. 4: The efficiency in dependence of E' = Ui^Mial L' 3 and R ut- The values of t\, t<i and ip have 
been optimized with help of the simplex method. Cycles violating the breakdown criteria have are 
plotted in red. 
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FIG. 5: Evaluation of the effect of material properties. The upper and lower rows show optimized 
efficiency r] and gain G, respectively, in dependency of relative permittivity e r and polymer con- 
ductivity o"dc- The graphs in the left column was for polymer modulus Y = 1.2 MPa, in the right 
for Y = 5 MPa. 
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FIG. 6: The efficiency r\ and the gain G in dependence of the Maxwell time rprj- The data are 
identical to Fig. [5] with Y = 1.2 MPa. The inset shows a semilogarithmic plot for r\. The big red 
diamonds show the results for the DEAP material. 



with e r . As has been suspected, but not clearly demonstrated before, it is found here that 
the stiffer material can achieve a higher efficiency, for identical ctdc an d s r . Gain G and 
efficiency 77 behave similarly, but now the stiffer material has the lower response, because 
they undergo smaller strains and thus smaller relative change in the capacitance. The drop 
in gain appears less pronounced than the increase in efficiency, which indicates that stiffer 
elastomers are advantageous for energy harvesting. 

Fig. shows also that both 77 and G depend roughly on e r x <j^q. We note that this 
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FIG. 7: Left panel: the efficiency w in dependence of the amplitude of the driving force /a and the 
stiffness Y = 8C10 of the material. Right panel: The efficiency and the gain in dependency of the 
A* = (JAmp/Y. The electrical material parameter are taken from[lj 

is nearly the expression for the time constant of the charge decay in the capacitor, rpc = 
e r 6o/a-Dc, also known as the Maxwell relaxation time. To elucidate this observation, rj and 
G are plotted against rpc, see Fig. |6j All efficiency values appear to fall below a particular 
threshold; below rpc = tc, the efficiency varies with rpc, above it is nearly constant and are 
distributed between 0.2 and 0.6. 

Also gain is seen to depend directly upon rpc- It collapses on a straight line in the 
semilogarithmic plot. Interestingly, it is only positive for rpc > tc, clearly showing that 
the Maxwell relaxation time must be shorter than the period of the mechanical driving. 
Comparing to the efficiency plot, high gain does not necessarily lead to high efficiency as 
was seen previously, see Fig. [3j In summary, the highest efficiencies and gains are slightly 
above 0.6 and encouragingly, the DEAP material shows both good efficiency and gain. 

These plots clearly show the importance of the proper choice of the time scales of the 
loss mechanisms, and that they must be chosen in accordance with the relevant driving time 
scale. Furthermore, rpc can be identified as the dominant material parameter for energy 
harvesting. 

Mechanical material parameter optimization 

The impact of the mechanical properties of the elastomer was evaluated by varying the 
stress-strain properties, simply by multiplication of Cio, C20, and C30 with the same factor. 
Then, the strength of the driving force was varied and either the efficiency rj or the gain 
G was optimized, as described above. As is clearly seen from Fig. [7] (left), the efficiency 
increases with higher stiffness, while it decreases for higher driving forces. 
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This data is efficiently analyzed when introducing a new variable, the virtual stretch 
A* = OAmp/y = fA/{L' 2 L' 3 /Y), which is a relative measure of the level of stretch during 
the cycle (see Fig. [7]). Remarkably, this plot shows that small stretches (originating from 
high stiffnesses or small driving forces) result in high efficiency. The gain behaves opposite, 
therefore energy harvesting can only be practicable in a narrow range of A*. 

CONCLUSION 

In this article, energy harvesting using electro active polymers has been studied by means 
of numerical simulations and analytical considerations. In detail, a realistic model describing 
the dynamics of the polymer as a relaxation process has been presented, as well as electrical 
control circuits and various limiting criteria. This system is driven by a periodic force, which 
compares well to an energy harvesting device using ocean waves or similar (fluid) mechanical 
forcing, as in generators in shoes or energy producing patches. The model is generic in that 
all its units describe essential and necessary parts of any future realization of an energy 
harvesting power plant. A direct coupling between the driving force and the soft capacitor 
was chosen, which applies well if complicated and expensive mechanical setups have to be 
avoided. 

The study focused on qualitative and quantitative measures (efficiency and gain) for 
energy harvesting and on the corresponding parameter optimization, which is important 
for the design of potential power plants. The study showed that optimizing charging and 
discharging times as well as the electrical loading parameters are crucial to obtain high 
efficiencies. Furthermore, it could be shown that the Maxwell time rpc and the virtual 
stretch A* are powerful tools in the analysis and should be tuned rationally to obtain optimal 
energy output. Among others, this leads to the strong conclusion that the Maxwell time 
must by larger than the period of the mechanical driving. 

The realistic simulation of the DEAP silicone material shows a specific harvesting en- 
ergy density of 5 Jkg -1 for a wave period of 10 s. Assuming a realistic raw material price 
for silicone of lOEuro/kg, this study indicates that a full return of materials investment is 
possible within 10 years. The investment will of course be compounded by other costs, for 
electronics, installations and maintenance. This time span is a very conservative estimate. 
It can be reduced if electrical fields larger than 20 V/zm -1 are used, which could be possible 



13 



with stiffer materials [32]. Further improvements should be possible with a more elaborate 
charging and discharging scheme that exploits the full variable range within the limit crite- 
ria \1'S\ ITT]. As such, this study shows that the wave energy harvesting technology based on 
soft capacitors has great potential for practical usability. 

APPENDIX: MECHANICAL MODEL 

The time dependent mechanical response of the polymer is given by Eq. (JT]) which de- 
scribes a Kelvin-like fluid. It assumes that shear stresses are not present and that the mate- 
rial will always relax back to the equilibrium. 7 is strongly related to the mechanical relax- 
ation time. In fact, the mechanical relaxation tm is given by r M x = dF(Ai)/dAi|A* ~ 
where F(Ai) is the RHS of @. A value 7 = 10 6 sPa" 1 was chosen resulting; in a relaxation 
time of tm ~ 1-8 s if external forces and charges are absent. 

To derive the final constitutive differential equation, the equations for Tn and T 33 are 
considered 

m M 1 E /l \ 

T u = a n + tr n - p = — - 7A1 

m M 1 E \ 

T33 = ^33 + ^33 - V = -7A3 , 

and the static pressure p is eliminated to obtain 

W - + K - <&) = ^ " ^ (Ai - A3) • (10) 

The force is applied in x-direction. Furthermore, T22 does not enter the game since pure- 
shear boundary conditions are assumed. 

The term af[ — crff = cxMech is the stress strain model and depends in general on the 
geometry of the polymer. The Yeoh model was chosen which agrees well with experimental 
results on the Danfoss material. It can be derived from an = XidW/dXi. W is the energy 
function which takes the form W = C w (h - 3) + C 20 (/i - 3) 2 + C 30 (h - 3) 3 for the Yeoh 
model with l\ = \\ + \\ + A|. Inserting the pure-shear assumptions and summarizing the 
terms for the x and the z directory lead finally to 

a M ech = 2 (\i - Ar 2 ) (C 10 + 2C 20 A + 3C 30 A 2 ) (11) 
with A = Af + 2 — 2. Usually the stiffness of materials is quantified by the Youngs modulus 
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E, which is defined as Y 
results in Y = 8Ci . 



OX U= 



. For the Yoeh model under pure shear conditions this 



The second term on the LHS of (10) is the Maxwell stress |24l [26] 

U 2 



C"Ele 



-e r £oE 3 — —e r eo—^\ 1 

L 3 



(12) 



where E 3 is the electrical field and U the voltage between the electrodes. 
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